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FOREWORD 


This  report  describes  results  from  computer  simulations  performed 
to  predict  the  effects  of  various  parameters  on  the  Beacon  Tracking 
System.  The  results  were  used  to  plan  the  implementation  of  the  Beacon 
Tracking  System  on  the  Naval  Surface  Weapons  Center/Dahl gren  Laboratory 
Test  Range.  Ted  Sims  of  DK-11  assisted  in  the  formulation  of  the 
computer  simulation  program.  This  report  has  been  reviewed  by 
R.  J.  Anderle,  Head,  Astronautics  & Geodesy  Division. 
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INTRODUCTION 


As  part  of  the  research  and  development  effort  in  producing  the 
hardware  and  software  for  the  Beacon  Tracking  System  , a simulation 
of  the  data  and  the  data  reduction  scheme  was  performed  to  investigate 
the  sensitivity  of  the  system  to  certain  key  parameters.  The  simulation 
consists  of  generating  synthetic  TDOA  data  for  a given  test  set  up, 
and  performing  a weighted  least  squares  solution  of  this  data  for  target 
position.  Target  maneuver  data  such  as  might  be  encountered  in  a true 
tracking  situation  is  included  in  the  simulation.  Each  of  the  features 
of  the  simulation  is  discussed  in  detail  below. 


SIMULATION 


1 . TDOA  Concept 

The  TDOA  concept  is  based  on  the  principle  that  a range  difference 
from  a spherical  wave  front,  centered  at  an  emitter,  to  two  known 
receivers  will  define  a hyperboloid  passing  through  the  emitter. 
Simultaneous  solution  of  three  hyperboloids  gives  three  dimensional 
coordinates  of  the  emitter.  Three  hyperboloids  are  defined  if  a 
minimum  of  four  receivers  (three  sets  of  two)  make  range  difference 
measurements.  Additional  sets  over-define  the  solution  permitting  in- 
vestigation of  additional  parameters. 

The  TDOA  system  studied  is  shown  in  Figure  1.  A transmitter  or 
beacon  emits  a pulse  that  travels  to  a master  receiver  and  several 
slave  transceivers  placed  at  known  locations.  Upon  receipt  of  the 
target  beacon  each  slave  transmits  a pulse  back  to  the  master.  Tm 
and  T]  are  the  times  the  target  pulse  travels  from  the  beacon  to  the 
master  and  a slave,  respecti vely.  Tim  is  the  time  a pulse  travels  from 
the  slave  to  the  master.  aT q is  the  slave  transceiver  reaction  time 
(hereafter  called  receiver  thru  time  or  RTT).  Then, 

TDOAi  = (T-j  + AT]  + Tim)  - Tm  (1) 

Since  the  locations  of  the  slave  and  master  stations  are  known,  T]m  is 
known.  It  is  the  distance  between  them  divided  by  the  speed  of  light. 
Receiver  thru  time  is  hardware  dependent  and  empirically  determined. 
Removing  the  two  known  quantities  leaves 

TDOA]  = T]  - TM  (2) 
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If  Xm>  Y^,  Zjv)  are  the  coordinates  of  the  master  receiver;  Xn,  Yn,  Zn 
are  the  coordinates  of  the  n slave  receivers,  and  X,  Y,  Z are  the  un- 
known beacon  position. 


tM  = £ [ (XM 

- x)2  + (YM  - Y)2  + (ZM  - Z)2]1/2 

(3) 

T,  - J-  [(X, 

- X)2  + (Y-]  - Y)2  + (Zt  - Z)2]1/2 

(4) 

Tn  = l [(Xn  - X)2  + (Yn  - Y)2  + (Zn  - Z)2]1/2  (5) 

For  (n+1)  receivers  (1  master,  n slaves),  n independent  TDOA  observations 
are  solved  simultaneously  for  X,  Y,  and  Z. 

The  n sets  of  simultaneous  equations  are  solved  by  the  method  of 
weighted  least  squares  discussed  in  Appendix  B. 


2.  Receiver  Geometry 

The  Beacon  Tracking  System  operates  in  the  geographic  area  shown 
on  the  map  in  Figure  2.  It  is  designed  to  track  over  the  range  from 
4 to  24  km.  The  master  station  is  located  at  the  main  firing  range  at 
NSWC/DL.  Slave  stations  are  located  along  the  shores  of  the  Potomac 
River,  down-range  from  NSWC/DL.  Three  criterion  are  met  in  the 
selection  of  each  of  the  slave  sites:  (1)  the  slave  station  has  a 

clear  line  of  sight  to  the  master  station  for  data  telemetry;  (2) 
the  geographical  survey  position  of  the  site  is  known  accurately,  and 
(3)  the  site  is  easily  accessable  to  NSWC/DL  personnel.  On  the  basis 
of  these  criterion  the  following  positions  are  chosen  for  study  in 
the  simulation:  station  5,  9,  25,  29,  12  (also  called  Swan  Point), 

Cobb  Island,  and  Oakland,  an  inland  point  on  the  Virginia  side  of  the 
river  (not  visible  in  Figure  2). 


3.  Coordinate  System 

For  the  Beacon  Tracking  simulation  the  various  slave  sites  are 
referenced  to  the  master  station  in  cartesian  coordinates  X,  Y,  Z. 

The  Master  station  is  at  0,0  with  positive  X "East",  positive  Y 
"North",  and  positive  Z "Up".  All  the  sites  are  surveyed  relative  to 
a reference  ellipsoid,  then  transformed  from  geodetic  to  cartesian 
coordinates  with  the  master  station  as  origin. 
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4.  Data 


Three  different  kinds  of  flight  paths  are  used  as  input  to  the 
tracking  simulation:  (1)  a straight  line  path  up  the  river  toward 
the  master  station  at  a velocity  of  about  450  mph  and  altitude  of  one 
kilometer  simulates  a nominal  aircraft  flight  path  (this  path  is 
called  a "riverflight"  for  reference);  (2)  an  aircraft  flight  path 
as  in  (1)  but  with  a 4 km  period  sine  curve  on  the  X,  Y,  and  Z com- 
ponents simulates  aircraft  maneuvers  (this  path  is  called  a 
"riversincurve" ) ; and  (3)  a straight  line  path  up  river  toward  the 
master  station,  velocity  about  250  mph,  with  an  abrupt  drop  in  altitude 
from  8000  ft.  to  4000  ft.  simulates  the  path  an  inbound  missile  might 
take.  This  path  is  called  a "missile  trajectory".  Using  these  three 
flight  paths,  synthetic  TD0A  data  is  generated  for  use  in  the  computer 
simulation. 

Sigma  on  the  data  class  is  taken  to  be  ao  = 1.X10-9.  A-priori 
estimate  of  the  hardware  (i.e.  receiver)  sigmas  is  taken  to  be 
°RTT  = 5.X10"9.  Atmospheric  refraction  effects  were  not  considered 
to  maintain  simplicity. 


RESULTS 


1 . Data  Noise 


To  test  the  effect  of  data  noise  on  the  least  squares  solution, 
a comparison  is  made  between  two  "riverfl ights" , identical  except  for 
noise  on  the  synthetic  data.  The  receiver  configuration  is  Master,  12, 
9,  25,  and  29.  Figure  3 shows  the  trace  of  the  sigma  levels  in  X,  Y, 
and  Z (sigmas)  over  the  path  from  32  km  to  2 km  with  no  noise  on  the 
data.  A second  "riverflight"  trajectory  is  run  identical  to  the  first 
except  that  0 to  +5  nanoseconds  of  random  noise  is  put  on  the  data  of 
the  second  trajectory.  Figure  4 shows  the  difference  between  the  true 
range  to  the  aircraft  and  the  range  computed  using  the  noisy  data. 

The  differences  in  range  shown  in  Figure  4 appear  significant  because 
the  solution  is  based  on  a single  data  burst  at  each  point.  The 
differences  can  be  greatly  reduced  by  obtaining  enough  data  to  permit 
an  "average"  of  several  bursts  of  data  to  be  taken  at  each  point. 

Note  that  in  both  Figures  3 and  4 an  abrupt  increase  in  solution 
error  is  seen  near  a range  of  26  km.  This  range  is  the  point  where 
the  aircraft  enters  an  unfavorable  geometry  spot  in  the  receiver 
configuration. 
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2.  Target  Maneuvers 

To  test  the  ability  of  the  Beacon  Tracking  System  to  follow  a 
maneuvering  target,  a "rivers incurve"  path  is  flown  using  the  receiver 
configuration  M,  12,  9,  25,  29.  ^nthetic  data  based  on  these 
factors  were  generated.  Figures  “and  6 show  the  effect  the  4 km 
period  sine  curves  on  X,  Y,  and  Z have-on  the  solution.  The  solution 
for  X and  Y (Figure  5)  show  no  wavering  over  t^fcarea  of  interest, 

4 to  24  km.  The  sine  path  is  seen  in  the  solution  for  Z (Figure  6); 
however,  acceptable  sigma  levels  are  maintained.  Figures  7,  8,  and  9 
show  the  displacement  in  X,  Y,  and  Z,  respectively,  between  the  "true" 
and  computed  positions  based  on  synthetic  data. 

The  reader  will  note  that  in  nearly  all  cases  presented  here, 
sigma  levels  in  the  solution  for  Z differ  from  X and  Y by  an  order 
of  magnitude.  This  difference  is  geometry  related  in  that  the  receiver 
stations  all  lie  in  approximately  the  same  plane.  This  "weakness" 
in  the  solution  for  Z would  not  be  seen  were  it  physically  possible 
to  position  the  receivers  in  three  dimensions  around  the  target  flight 
path. 

Figure  10  shows  results  from  a case  that  is  a duplicate  of  the 
above  "riversincurve"  case  (Figures  5 and  6)  with  one  small  change. 

In  Figure  10  the  least  squares  solution  is  computed  using  five  iterations 
instead  of  three  as  in  the  previous  case.  No  effect  is  demonstrated 
from  the  five  iteration  case  over  three  iterations.  Thus, all  other 
cases  are  solved  using  three  iterations. 

To  test  the  ability  of  the  Beacon  Tracking  System  to  track  the 
path  of  an  incoming  missile,  a "missile  trajectory"  is  flown  through 
the  configuration  M,  12,  9,  25,  29.  The  target  travels  at  250  mp'ri  at 
an  altitude  of  8,000  ft.  until  it  reaches  a range  of  13  km.  Then  it 
rapidly  drops  to  an  altitude  of  4,000  ft.  Synthetic  data  is  generated 
from  these  conditions.  Figure  11  shows  no  drastic  degradation  in  the 
solution  for  X and  Y.  Figure  12  shows  a slight  degradation  in  the 
solution  for  Z at  the  altitude  drop,  but  with  quick  recovery.  Figures 
13,  14,  and  15  show  the  displacement  between  the  "true"  and  the 
computed  values  for  target  position  X,  Y,  and  Z,  respectively.  No 
loss  of  accuracy  is  seen  in  X and  Y;  barely  perceptible  accuracy 
changes  are  seen  in  Z. 


3.  Receiver  Positions 

* Our  criterion  that  the  receiver  sites  be  easily  accessible  to 
NSWC/DL  personnel  make  it  desirable  to  have  slave  sites  on  the  Virginia 
side  of  the  Potomac  River  only.  To  test  the  effect  on  the  system 
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solution  if  no  slave  receivers  are  used  on  the  Maryland  side,  two  cases 
are  given  below  for  two  different  configurations. 

In  the  first  case  where  no  Maryland  receivers  were  used  a "river- 
sincurve"  is  flown  using  receiver  configuration  M,  9,  25,  29.  Synthetic 
data  were  generated  with  these  conditions.  Figures  16  thru  20  show  the 
solution  sigma  levels  and  the  displacements  in  X,  Y,  and  Z between 
"true"  and  computed  positions.  The  solution  for  X in  this  configuration 
is  acceptable  over  the  range  of  interest,  4 to  24  km.  The  solutions 
for  Y and  Z,  however,  are  not.  Both  Y and  Z show  poor  accuracy  level 
over  the  range  from  16  to  24  km.  In  addition  Y and  Z solutions  degrade 
significantly  around  23  km.  This  degradation  is  attributed  to  the 
receiver  geometry.  Because  of  poor  solution  confidence  in  two  dimensions 
and  the  solution  degrading  within  our  area  of  interest,  this  configuration 
is  judged  inadequate  for  Beacon  Tracking. 


The  second  case  where  no  Maryland  receivers  were  used  is  a "missile 
trajectory"  with  receiver  configuration  M,  5,  9,  25.  Synthetic  data 
were  generated  with  these  conditions.  Figures  21  thru  25  show  the 
solution  confidence  levels  and  the  X,  Y,  Z displacements  between  "true" 
and  computed  target  positions.  Substituting  station  5 for  station  29 
eliminates  the  solution  degradation  due  to  receiver  geometry  within  our 
range  of  interest.  The  solution  is  acceptable  in  X and  Y for  this 
configuration,  but  questionable  in  Z. 


We  add  two  receivers  in  Maryland  to  the  above  case  to  test  solution 
improvement.  The  configuration  is  now  M,  5,  9,  25,  12,  Cobb.  Figures 
26  thru  30  show  the  new  solution  based  on  data  from  this  new  configura- 
tion with  the  Maryland  sites  added.  Over  the  range  of  interest  solution 
improvement  is  dramatic.  In  X and  Y solution  accuracy  improves  by  a 
factor  of  two  or  more.  In  Z,  however,  solution  accuracy  improves  by  a 
factor  of  four  or  more.  Improvement  in  Z is  particularly  seen  in  the 
range  from  14  km  to  24  km.  Where  the  Z solution  is  questionable 
in  this  range  for  the  no  Maryland  receivers  case,  it  is  clearly 
acceptable  with  the  Maryland  receivers  added. 


It  is  recognized  that  two  factors  affect  the  solution  improvement 
shown  in  the  above  case:  (1)  the  addition  of  two  more  receivers  to 

the  configuration,  and  (2)  the  positioning  of  receivers  "around"  the 
target.  Other  NSWC  analysis  has  demonstrated  that  the  second  factor 
is  the  most  significant.  This  configuration,  M,  5,  9,  25,  12,  Cobb 
is  the  one  selected  for  the  operational  Beacon  Tracking  System. 
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4.  Receiver  Failure 


At  the  time  of  this  simulation  study  the  Beacon  Tracking  System 
hardware  was  experimental  and  its  reliability  as  yet  unproven. 
Therefore,  the  possibility  existed  that  one  or  more  slave  receivers 
would  not  function  during  the  course  of  a tracking  sequence.  To  test 
the  sensitivity  of  the  system  solution  to  the  loss  of  a receiver,  two 
such  scenarios  are  simulated. 

In  the  first  a "riversincurve"  is  flown  with  receiver  net  M,9,25, 

12,  Cobb,  Oakland.  At  each  data  point  each  receiver  is  alternately 
dropped  and  the  solution  performed  without  that  receiver's  data  con- 
tributing to  the  solution.  Figures  31  and  32  show  the  resultant  plots 
of  sigma  levels  for  X,Y  and  Z,  respectively.  For  this  configuration 
stations  9 and  25  are  shown  to  be  the  most  critical  to  the  solution 
confidence.  Station  9 is  particularly  critical  in  all  three  dimensions. 
Without  data  from  Station  9,  the  solution  degrades  drastically  in 
the  range  from  12  km  to  24  km.  Station  25  is  critical  to  the  solution 

for  X and  Y in  the  distant  ranges  (24+  km).  The  other  three  stations, 

Oakland,  Swan,  and  Cobb,  are  of  about  equal  but  lesser  importance  in 
X and  Y.  Station  9 is  the  critical  one  in  the  solution  for  Z,  with  25, 

Swan,  Oakland,  and  Cobb  of  about  equal  but  far  lesser  importance. 

In  the  second  scenario  a "riversincurve"  is  flown  with  receiver 
configuration  M,  5,  9,  25,  12,  Cobb.  As  the  geometry  suggests,  station 
25  is  again  the  most  critical  to  the  solution  in  all  three  dimensions, 
particularly  at  the  distant  ranges.  For  X and  Y (see  Figure  33)  the 
other  four  slave  stations  are  of  about  equal  but  far  lesser  importance. 
Of  particular  interest  for  this  scenario  is  the  importance  of  the 
various  slaves  in  the  solution  for  Z,  as  shown  in  Figure  34.  Again, 
station  25  is  the  most  critical  to  the  solution.  However,  stations  9 
and  12  (Swan)  are  significantly  more  important  to  the  Z solution  than 
are  stations  5 and  Cobb.  This  reconfirms  our  earlier  observations 
about  the  importance  of  at  least  one  slave  site  on  the  Maryland  side 
of  the  river. 
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Z DISPLACEMENT  VS.  RANGE,  MISSILE  PATH 
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X AND  Y SIGMA  VS.  RANGE,  SINE  PATH,  NO  MARYLAND 
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Z SIGMA  VS.  RANGE,  SINE  PATH,  NO  MARYLAND 
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X DISPLACEMENT  VS.  RANGE,  SINE  PATH,  NO  MARYLAND 
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Y DISPLACEMENT  VS.  RANGE,  SINE  PATH,  NO  MARYLAND 
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Z DISPLACEMENT  VS.  RANGE,  SINE  PATH,  NO  MARYLAND 
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Z SIGMA  VS.  RANGE,  MISSILE  PATH,  NO  MARYLAND 
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DISPLACEMENT  VS.  RANGE,  MISSILE  PATH,  NO  MARYLAND 
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Y DISPLACEMENT  VS.  RANGE,  MISSILE  PATH,  NO  MARYLAND 
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Assume,  as  in  the  discussion  on  page  2,  we  are  able  to  resolve  a 
quantity  T.  Let  these  measured  values  of  T be  denoted  by  TOBS.  Assume 
the  quantity  measured  is  a function  of  some  parameter  space  T = T( Pi , . . . » ) . 
Of  interest  are  the  specific  values  of  the  parameter  set  P-|,...,Pk  that 
relate  to  the  measured  values  TOBS. 

Taking  M measurements  of  T results  in  the  following  system  of 
equations,  relating  measured  values  of  T to  the  parameters  P-|,...,Pk 

(TOBS ) -j  = T-|(Pi Pk)  for  i = 1,2,...,M  (1) 

If  M>K,  the  system  is  overdetermined,  i.e.  there  are  more  equations  than 
unknown  parameters.  In  order  to  obtain  a unique  solution  to  Equation 
1,  an  additional  criterion,  the  least  square  condition,  is  imposed. 

This  criterion  is  that 


3 [ (TOBS ) i ) - Ti (Pi Pk)]2 


be  a minimum,  ai  is  the  standard  deviation  of  the  ith  measurement.  Let 
(P-|0 Pk0)  be  an  initial  estimate  to  the  desired  (Pi , . . . Pk ) • Let 

(TCOMP)i  = Ti(Pi Pk)  (2) 

Expanding  T(Pi , . . . ,Pk)  in  a Taylor  Series  around  the  point  (P]0,. . . ,Pk°) 
results  in: 

k 3T  ■ 

Ti(P1,...,Pk)  = T-j  (P]  °, . . . ,Pk°)  + Z — - (Pj-Pj°)  + higher  terms  (3) 

j = l 3\j 

by  neglecting  the  higher  order  terms  and  defining 

APj  = (Pj  - Pj°)  (4) 


the  least  squares  condition  is  satisfied  if  for  all  parameters 
PL(L  = 1.....K) 


9 M 
9P 


[(TOBS)-j  - (T-j  (P-,^,. . . ,Pk  ) ~ 


ti  3P.i 


ill  APj)]2 


i = l 


= 0 
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Let  (0-C) i = (TOBS) i - ^ (P1,...,Pk).  Then 


M (O-C)j2 
i = l a 9PL 


k M i »T-  3Ti 

Z Z — • - — — AP, 

j=l  i = l oi2  3Pj  3PL  J 


(6) 


There  is  one  such  expression  for  each  parameter  Pi  . Once  the  initial 
guess  (PL0,...,Pk0)  is  given,  all  quantities  in  tne  above  expression 
are  known  except  for  APj . Solving  for  APj  results  in  what  is  termed 
an  improvement  to  the  jth  parameter. 

The  above  equations  can  be  written  in  matrix  form  as  follows: 

. = 1 2li 

Aik  " qT  8 P k 

1 3^ i 3^ i 

Bik  j ^i7  ^k 

Note  that  B is  symmetric  and  that  B = A^A.  Equation  6 becomes 

(O-C)i2 

? ^ Aik  T:  Bkj  APj- 

l J 
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I 
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Let 


Aik* 


/APl\ 

and  aP  =f  : \ . 

\ / 

Then  BaP  = E (7) 

Solution  of  the  normal  equations  (7)  results  in  AP  or  the  improvements 
to  the  parameter  set  (P-|0,. . .Pk0).  if  the  initial  estimate  (P-|0, . . . ,PkB) 
is  sufficiently  good,  defining 

(P! ,Pk)  = (P-,0  + A?! pk°  + APk) 

results  in  the  best  (in  the  least  squares  sense)  solution  to  the  minimi- 
zation of 


B-l 


N [(TOBS)i  -Ti  (Pls...,Pn)]2 


Because  we  have  made  a linear  approximation  to  T(Pi , . . . ,P|<) , we  have 
minimized  a linear  problem.  Since  in  general  T - T(P-| ,. . . ,P|<)  is  not 
linear,  it  is  necessary  to  iterate  on  this  procedure  until  a best 
solution  is  determined. 

In  order  to  achieve  an  over-determined  system,  we  make  use  of 
a-priori  knowledge  of  the  standard  deviation  of  some  of  the  parameters. 
The  uncertainties  in  the  RTT  are  determined  by  independent  calibrations. 
Thus,  after  we  form  the  system 


BaP  = E, 


to  B we  add  the  diagonal  matrix 


where  ai  is  the  standard  deviation  of  the  ith  parameter.  The  above 
diagonal  matrix  contains  n non-zero  terms,  ox,  ay,  oz  (the  standard 
deviation  of  the  position  of  the  beacon  transmitter)  is  assigned  an 
infinite  "sigma"  and 


1 and  and  _J_  are  all  zero. 

? 2 2 
ox^  0/  oz‘ 


Now  we  have  (with  n receivers)  3 + n parameters,  n - 1 observations,  and 
n a-priori  observations.  For  r >4,  n-l+n>3+n  and  the  system  is  over- 
determined. 

The  above  diagonal  matrix  (i.e.  the  matrix  of  a-priori  estimates 


B-2 


on  RTT's)  is  added  to  the  left-hand  side  of  the  normal  equations 
only  once.  This  one  time  occurs  at  the  first  point  at  which  the 
normal  equations  are  formed. 
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